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Abstract. The origin of the deficiency in neutral Hydrogen of 13 spiral galaxies lying in the outskirts of the Virgo 
cluster is reassessed. If these galaxies have passed through the core of the cluster, their interstellar gas should have 
been lost through ram pressure stripping by the hot X-ray emitting gas of the cluster. We analyze the positions 
of these Hi-deficient and other spiral galaxies in velocity-distance plots, in which we include our compilation of 
velocity-distance data on 61 elliptical galaxies, and compare with simulated velocity-distance diagrams obtained 
• from cosmological iV-body simulations. We find that ~ 20% relative Tully-Fisher distance errors are consistent 

with the great majority of the spirals, except for a small number of objects, whose positions in the velocity-distance 
diagram suggest grossly incorrect distances, which implies that the Tully-Fisher error distribution function has 
non-gaussian wings. Moreover, we find that the distance errors may lead to an incorrect fitting of the Tolman- 
Bondi solution that can generate significant errors in the distance and especially the mass estimates of the cluster. 
We suggest 4 possibilities for the outlying Hi-deficient spirals (in decreasing frequency): 1) they have large relative 
distance errors and are in fact close enough (at distances between 12.7 and 20.9 Mpc from us) to the cluster to 
have passed through its core and seen their gas removed by ram pressure stripping; 2) their gas is converted to 
stars by tidal interactions with other galaxies; 3) their gas is heated during recent mergers with smaller galaxies; 
and 4) they are in reality not Hi-deficient (e.g. SO/a's misclassified as Sa's). 

Key words. Galaxies: clusters: individual: Virgo - galaxies: evolution - galaxies: ISM - cosmology: distance scale 
- methods: N-body simulations 



1. Introduction 

During the last thirty years, substantial work has been 
carried out to quantify the effects of the environment on 
the properties of galaxies, i.e. to distinguish the character- 
istics of galaxies in dense, aggressive regions such as clus- 
ters, especially cluster cores, and low density zones away 
from clusters and groups. The nearby Virgo cluster is an 
excellent candidate for such studies, because its proxim- 
ity allows 1) an easy comparison with field galaxies and 
2) the measurement of distances independent of radial ve- 
locity, providing three-dimensional information on cluster 



properties (e.g.lYoung fc Curriejl995tlGavazzi et al . 1999; 



ISolanes et ail 2002) impossible to obtain in other clusters. 



One particular characteristic of spiral galaxies is their 
neutral Hydrogen (HI) gas content. It has been known 
for some time that a class of H I- deficient galaxies ex- 
ists in the centers of clusters, such as t he Virgo cluster 
( Chamaraux. Balkowski. fc Gerardlll980|) . and their defi- 
ciency, normalized to their optical diameter and morpho- 
logical tv_p_e J Js_anti^orre]ated_v^ dis- 
tance (|Havnes fc Giovanellilll98ri ISolanes et alJl200l[) . 
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IChamaraux et al.1 l)l980|) considered ram pressure 
stripping of interstellar atomic Hydroge n by the hot 
diffuse intracluster gas, first proposed by iGunn fc Gotil 
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as one cause of the H I deficiency pattern of Virgo 
spirals. For galaxies moving face-on, ram pressure scales 
as the intracluster gas density times the square of the 
galaxy's velocity relative to this intracluster gas, so the 
role of ram pressure stripping was strengthened by the 
correlation of HI deficiency with the X-ray luminosity of 
clusters found bylGiovanelli &: Havnesl ( 19851 see however 
ISolanes et alJl200 1) and by the presen ce of truncated H I 
disks observed by Cavatte et al.1 1 19901) . 

It hence came as a surprise when lSolanes et al 
discovered spiral galaxies, whose Hydrogen content was 
deficient relative to other spiral galaxies of the same mor- 
phological type, far away from the core of Virgo, with 
several deficient spirals in the foreground of Virgo and 
particularly one group of spirals 7 to 10 Mpc (depend- 
ing on the adopted distance to Virgo) behind the cluster 
with one-third of its 15 spirals displaying HI deficiencies 
deviating by more than 3 a from normalcy. 

In a followup study, ISanchis et ail l|2002l) explored the 
possibility that the group of deficient spirals behind Virgo 
might have plunged into the cluster core a few Gyr ago 
and lost their gas by interacting with the hot intraclus- 
ter medium. Sanchis et al. applied a simple dynamical 
model that considered not only the classical deduction of 
the virgoccntric infalling regime, but also the trajectories 
of galaxies that had already crossed once the center of 
the core. This model succeeded in explaining the overall 
characteristics of the observed velocity-distance diagram, 
including the envelope of the virialized cluster drawn by 
galaxies on first rebound. With the available data, they 
argued, from a statistical point of view, that one could 
not discard the possibility that some of the gas-deficient 
objects, including the background group in the outskirts of 
the Virgo Cluster Region, were not newcomers, but have 
passed through the cluster and lost their gas by ram pres- 
sure str ipping. 

We ijMamon et alJ Ep04) have recently estimated the 
maximum radius out to where galaxies falling into a clus- 
ter can bounce out. We found that this maximum rebound 
radius corresponds to 1 to 2.5 cluster virial radii. After es- 
timating the virial radius of the Virgo cluster using X-ray 
observations, we concluded that the Hi-deficient galaxies 
located in the foreground or background of the Virgo clus- 
ter appear to lie much further from the center of the clus- 
ter than these maximum distances, and therefore that ram 
pressure stripping cannot account for the H I-dcficicncy of 
these outlying galaxies if their distances are correct. 

In this paper, we investigate in further detail the ori- 
gin of the H I-deficiency in the spiral galaxies that appear 
too far from the Virgo cluster to have suffered ram pres- 
sure stripping of their interstellar gas, by 1) adding to 
our four-dimensional (angular position, distance and ra- 
dial velocity) sample of spiral galaxies, our compilation of 
the literature for elliptical galaxies towards Virgo, leading 
us to identify more precisely the density of the environ- 
ment of the background deficient galaxies and estimate 
more accurately the distance to the center of Virgo, which 
is a crucial ingredient for the dynamical infall model, and 



2) using cosmological iV-body simulations that have the 
advantages of their lack of observational errors, complete- 
ness, and of defining a realistic infall pattern. 

The sample of elliptical galaxies is described in Sect.FJl 
In Sect. |31 after describing the cosmological TV-body sim- 
ulations used here, we compare the simulated and obser- 
vational samples to evaluate the influence of errors in dis- 
tance estimates. Finally in Sect.|U we discuss possible ex- 
planations for the HI deficiency of the outlying galaxies. 

2. The inclusion of ellipticals into the Virgo 
galaxy sample 

Using only late- type spirals. ISanchis et alJ l|2002j) applied 
a spherical infall model to the velocity field of the Virgo 
Cluster Region and obtained two possible distances of 
20 and 21 Mpc. Other distance estimates to the Virgo 
cluster, based upon spir al galaxies, produced distances in 
the range from 16 M pc l|Tullv fc Shavalll984h to 21 Mpc 
llEkholm et al.ll200(t) . 

An alternative way to estimate the cluster distance 
is by using the early-type components of the cluster, as 
they are thought to be better tracers of the center of 
the cluster due to morph ological segregation (see e.g. 
iFerguson fc Binggelil ll994). To do so, we have collected 
as many as possible early-type galaxies (E's, dE's and 
SO's) with redshift-independent distance measurements. 
Heliocentric velocities were extracted from the LEDA cata- 
log of the HyperLEDA database and converted to the Local 
Group reference frame, while the direct distant measure- 
ments are obtained from a variety of sources. 

Table 1. Datasets contributing to elliptical sample 



Method 


Authors" 


Band 


rel. accuracy 


Number 








(rms) 




FP 


G 


H 


0.41 


41 


L — n 


Y L 


Bj 


0.54 


9 


R-n 


Yr 


Bj 


0.54 


14 


SBF 


T 


I 


0.18 


31 


SBF 


h 


I 


0.07 


6 


SBF 


Jif 


K' 


0.17 


7 


SBF 


N 


I 


0.15 


15 


GCLF 


K 


V-I 


0.07 


8 



a FP = Fundamental Plane; L — n, R — n — Luminosity / 
Sersic shape and Radius / Sersic shape relations; SBF = 
Surface Brightness Fluctuations : GCLF = Gl o bular Cluster 
Luminosity F unctio n. b G = iGavazzi et alJ (I 1999T) : T = 
iTonrv et, all l)200ll); Ji,Jk =~ l.lensen. Tonrv fc Lunninol 

1 20001): K 



1119981): N = [ Neilsen fc Tsvetangvl 
iKundu fc Whitmorel (1200 ll) : Yz,,Y H = 
(1993). 



Young fc Currid 



Tablc^summarizes the distance datasets on the early- 
type galaxies towards Virgo. Although different authors 
use different methods with different accuracy, we have 
not found systematic deviations between them, and fur- 
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thcrmore found that the spread in differences of the dis- 
tance moduli was consistent with the quoted errors, for 
all sets of early- type galaxies common to two authors. 
Therefore, we estimated a global distance modulus d for 
each galaxy, using the maxim um likelihood estimator (e.g. 
iBevineton fc RobinsorJll992() 



d = 



(i) 



where di are the distance moduli and a% are the abso- 
lute errors on these distance moduli. Equation (JJJ sup- 
poses that the relative error on distance measurements 
is independent of distance and that the p.d.f. of the 
distance moduli (proportional to log distance) is gaus- 
sian. The error on the global distance modulus is (e.g. 
iBevington fc R.obinsorJ) 



a4 = 



(2) 



Although this estimator of the error on the distance mod- 
ulus is independent on the scatter of the measurements, 
we show in Appendix 1X1 (eas. |A.7j and |A.8j ) that equa- 
tions (JIJ and J5J also represent the median and the half- 
width (as defined by (x$4 — Xie)/2) of the likelihood func- 
tion, respectively. 

Table [21 provides our final sample of 63 early-type 
galaxies. Columns (1), (2) and (3) give the name of the 
galaxy, with its coordinates in columns (4) and (5). In col- 
umn (6), we give the distance (in Mpc, from eq. [Q) and 
its relative error (eq. 0) in column (7). Column (8) gives 
the velocity relative to the Local Group and column (9) 
provides the sources for the distance. 

Figure ^ shows a wedge diagram of the Virgo region 
with the 61 early-types shown as circles and the 146 spirals 
shown as triangles. One clearly notices the concentration 
of elliptical galaxies at a distance of D ~ 17 Mpc, which 
is spread out in velocity space. 

3. Observed and simulated velocity-distance 
relation towards the Virgo Cluster 

For a better understanding of the kinematics of galaxies 
in the direction of the Virgo cluster, we have simulated an 
observation of the velocity-distance relation using the dark 
matter particles in the GallCS cosmological simulation. 

3.1. N-body simulations 

The iV-body simulations used here to fin d TV-body rep lica 
of the Virgo cluster were carried out bv iNininl l)l999|) . as 
described by lHatton et alJ l|2003j) in the context of the 
GallCS hybrid AT-body/semi-analytic model of hierarchi- 
cal galaxy formation. Here, we are only interested in the 
density and velocity fields directly traced by dark matter 
particles. The TV-body simulation contains 256 3 particles 
of mass 8.3 x 10 9 M Q in a box of 150 Mpc size and it is 
run with a softening length amounting to a spatial resolu- 
tion of 29 kpc. The simulation was run for a flat universe 



with cosmological parameters f2o = 0.333, Q\ = 0.667, 
Hq = 66.7kms~ 1 Mpc -1 , and erg = 0.88. Once the sim- 
ulation is run, halos of dark matter are detected with a 
'Friends-of-Friends' (FoF) algorithm 



flDayjS fit al -1 11? 

with a variable linking length such that the minimum mass 
of the FoF groups is 1.65 x 1O 11 M (20 particles) at any 
time step. With this method, over 2 x 10 4 halos are de- 
tected at the final timestep, corresponding to the present- 
day (z — 0) Universe. The GallCS halo finder does not 
allow halos within halos, so that a cluster, to which is as- 
signed a massive halo, cannot contain smaller halos within 
it. 



3.2. Picking isolated halos 

Among the 12 halos that GallCS produces at z — with 
M FoF > 10 14 M Q , we wish to choose those that resemble 
the Virgo cluster, both in terms of mass and environment. 

As we can see in the radial phase space diagrams 
of dark matter ha los of the simulations (Fig. 1 of 
IM anion et al.l l2004[) , it is quite common to find small 
groups around the main cluster halo, whereas the Virgo 
cluster is not known to have massive neighbors, except for 
the M49 group just within the virial radius. Note that the 
size of the simulation box, L — 150 Mpc, implies, given 
the periodic boundary condition, a maximum separation 
between halos of V^(L/2) = 130 Mpc. 

The effect of a neighboring halo will be from its tide 
on the test halo. The ratio of the tidal acceleration on a 
particle of the test halo at distance r from its center caused 
by a neighboring halo at distance R from the center of 
the first halo to the acceleration of this particle by the 
potential of its halo amounts to a mean density criterion: 



GM 2 {R) 
atid ~ — r 



and 



fthalo — 



R ? > 



GMi(r) 



so that 



(hid 
Ohalo 



P2(R) 

pi(r) 



(3) 



(4) 



(5) 



where the densities are average (and not local) and sub- 
scripts 1 and 2 stand for the central and the neighboring 
halo respectively. For clarity, we denote r 1; Mi, r 2 and 
Mi the virial radii and masses for the test and neighbor- 
ing halo, respectively (i.e. we drop the TOO' subscript). 
Then, equating r = n, noting that the mean density at 
the virial radius is the same for each halo, i.e. 100 times 
the critical density of the Universe, and writing 



R 

n= — 

T2 



(6) 



the influence of each neighbor is measured by the mean 
density of the neighbor measured at the center of the test 
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Table 2. Final sample of early- type galaxies 





Galaxy 




RA 


Dec 


D 


rel. error 


V 


sources 


vec 


NGC 


Messier 


(J2000) 


(Mpc) 




(kms- 1 ) 




49 


4168 




12 h 12 m 17!3 


+13°12'17" 


34.20 


0.41 


2128 


G 


220 


4233 


— 


12 h 17 m 07!7 


+07°37'26" 


32.96 


0.41 


2181 


G 


312 


4255 


— 


12 h 18 m 56"l 


+04°47'09" 


28.71 


0.41 


1603 


G 


319 


— 


— 


12 h 19 m 02!0 


+13°58'48" 


7.59 


0.38 


-266 


Y L Y R 


342 


4259 


— 


12 h 19 m 22!2 


+05°22'34" 


57.02 


0.41 


2306 


G 


345 


4261 


— 


12 h 19 m 23!2 


+05°49'32" 


32.27 


0.16 


2012 


GT 


355 


4262 


— 


12 h 19 m 3Cfi6 


+14°52'38" 


14.86 


0.41 


1181 


G 


369 


4267 


— 


12 h 19 m 45!4 


+12°47'53" 


14.66 


0.41 


852 


G 


575 


4318 


— 


12 h 22 m 43?!3 


+08°11'52" 


25.35 


0.41 


1047 


G 


648 


4339 


— 


12 h 23 m 34l9 


+06°04'54" 


16.44 


0.17 


1108 


T 


685 


4350 


— 


12 h 23 m 57!7 


+16°41'33" 


14.86 


0.41 


1049 


G 


731 


4365 


— 


12 h 24 m 28f!3 


+07° 19' 04" 


22.07 


0.05 


1058 


J/JkNT 


763 


4374 


84 


12 h 25 m 03!7 


+12°53'13" 


17.71 


0.07 


733 


GNT 


778 


4377 


— 


12 h 25 m 12!3 


+14°45'43" 


15.49 


0.41 


1182 


G 


784 


4379 


— 


12 h 25 m 14'8 


+15°36'26" 


14.19 


0.38 


871 


T 


810 


— 


— 


12 h 25 m 33"8 


+13°13'30" 


12.36 


0.38 


-516 




828 


4387 


— 


12 h 25 m 41!8 


+12°48'35" 


17.21 


0.35 


373 


GT 


881 


4406 


86 


12 h 26 m 12!2 


+12°56'44" 


17.40 


0.04 


-465 


J/JkKNT 


940 


— 


— 


12 h 26 m 47:i 


+ 12°27'14" 


7.24 


0.54 


1337 


Yr 


953 


— 


— 


12 h 26 m 54 ! :7 


+13°33'57" 


18.88 


0.38 


-678 


Y l Y r 


965 


— 


— 


12 h 27 m 03 ! !l 


+12°33'39" 


9.29 


0.54 


614 


Yr 


1025 


4434 


— 


12 h 27 m 3tf7 


+08°09'14" 


25.50 


0.15 


897 


GT 


1030 


4435 


— 


12 h 27 m 40 f :6 


+13°04'44" 


13.12 


0.41 


611 


G 


1146 


4458 


— 


12 h 28 m 57!6 


+13°14'29" 


17.97 


0.09 


499 


GKNT 


1173 


— 


— 


12 h 29 m 14!8 


+12°58'41" 


14.96 


0.38 


2292 


Y L Y R 


1196 


4468 


— 


12 h 29 m 30:9 


+14°02'55" 


15.85 


0.13 


732 


T 


1226 


4472 


49 


12 h 29 m 46"7 


+07°59'59" 


16.27 


0.03 


693 


GJ/JkKNT 


1231 


4473 


— 


12 h 29 m 48?!9 


+13°25'46" 


17.12 


0.05 


2041 


GKNT 


1242 


4474 


— 


12 h 29 m 53"7 


+14°04'06" 


15.14 


0.41 


1405 


G 


1250 


4476 


— 


12 h 29 m 59"l 


+12°20'55" 


19.08 


0.11 


1767 


NT 


1253 


4477 


— 


12 h 30 m 02:2 


+13°38'11" 


16.44 


0.41 


1164 


G 


1279 


4478 


— 


12 h 30 m 17!5 


+12°19'40" 


16.68 


0.09 


1220 


GNT 


1297 


4486B 


— 


12 h 30 m 32!0 


+12°29'27" 


16.31 


0.03 


136 


KN 


1308 


— 


— 


12 h 30 m 4ff!9 


+11°20'36" 


17.86 


0.38 


1511 


Y L Y R 


1316 


4486 


87 


12 h 30 m 49"4 


+ 12°23'28" 


16.82 


0.09 


1095 


GNT 


1321 


4489 


— 


12 h 30 m 52 ! !3 


+16°45'30" 


17.04 


0.13 


772 


JxT 


1386 


— 


— 


12 h 31 m 51"5 


+12°39'24" 


7.08 


0.38 


1135 


Y L Y R 


1412 


4503 


— 


12 h 32 m 06°3 


+iri0'35" 


11.75 


0.41 


1184 


G 


1420 


— 


— 


12 h 32 m 12"3 


+12°03'42" 


11.75 


0.38 


847 


Y l Yr 


1489 


— 


— 


12 h 33 m 14'0 


+10°55'43" 


11.30 


0.38 


-94 


Y l Yr 


1535 


4526 


— 


12 h 34 m 03"l 


+07°41'57" 


16.61 


0.17 


427 


GT 


1537 


4528 


— 


12 h 34 m 06!2 


+11°19'15" 


14.52 


0.41 


1195 


G 


1539 


— 


— 


12 h 34 m 06f3 


+12°44'40" 


8.83 


0.54 


1216 


Yr 


1549 


— 


— 


12 h 34 m 1419 


+ir04'16" 


19.32 


0.54 


1202 


Yr 


1619 


4550 


— 


12 h 35 m 3Cfi7 


+12°13'13" 


17.84 


0.09 


208 


KNT 


1630 


4551 


— 


12 h 35 m 38°l 


+12°15'49" 


17.24 


0.15 


973 


GT 


1632 


4552 


89 


12 h 35 m 40!0 


+ 12°33'22" 


15.31 


0.04 


115 


GJ/JkKNT 


1664 


4564 


— 


12 h 36 m 27!0 


+11°26'18" 


14.86 


0.15 


947 


GT 


1669 




— 


12 h 36 m 30:6 


+13°38'18" 


8.79 


0.54 


427 


Yr 


1692 


4570 


— 


12 h 36 m 53!5 


+07°14'46" 


13.80 


0.41 


1559 


G 


1720 


4578 




12 h 37 m 30:6 


+09°33'16" 


17.94 


0.06 


2102 


GJ/JkT 


1834 


4600 




12 h 40 m 23!0 


+03°07'02" 


7.35 


0.20 


615 


T 


1869 


4608 




12 h 41 m 13M 


+10°09'18" 


20.23 


0.41 


1633 


G 


1883 


4612 




12 h 41 m 32!8 


+07°18'50" 


14.06 


0.41 


1697 


G 


1902 


4620 




12 h 41 m 59!4 


+12°56'32" 


21.28 


0.28 


1012 


T 


1903 


4621 


59 


12 h 42 m 02!3 


+11°38'45" 


14.80 


0.05 


258 


GKNT 
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Galaxy 




RA 


Dec 


D 


rel. error 


V 


sources 


vcc 


NGC 


Messier 


(J2000) 


(Mpc) 




(kms -1 ) 




1938 


4638 


— 


12 h 42 m 47!5 


+ 11°26'32" 


17.33 


0.21 


956 


GT 


1939 


4624 




12 h 42 m 5010 


+02°41'16" 


15.14 


0.06 


924 




1978 


4649 


60 


12 h 43 m 40 s 1 


+11°33'08" 


16.36 


0.04 


969 


GKNT 


2000 


4660 




12 h 44 m 32 s 


+ 11°11'24" 


15.66 


0.11 


914 


GNT 


2087 


4733 




12 h 51 m 06":8 


+ 10°54'43" 


14.93 


0.18 


772 


T 


2092 


4754 




12 h 52 m 1715 


+11°18'50" 


16.94 


0.12 


1206 


GT 


2095 


4762 




12 h 52 m 55!9 


+11°13'49" 


10.81 


0.41 


818 


G 




Fig. 1. Positions of 146 spiral (triangles) and 61 elliptical (circles) galaxies in the Virgo field (within 9° from M87) 
as a function of Declination and distance (assumed correct, top) and radial velocity (bottom). Hi-deficient galaxies 
(deviating more than 3<r from normalcy) are shown as filled triangles — with the ones apparently outside the cluster 
(Table 7, below) shown as surrounded filled triangles (red triangles in the Electronic version of the Journal), while 
Messier ellipticals (see Table 01 are shown as filled circles, with M49 to the South, and M59, M89, M60, M87, M86 
and M84 at increasing distances (top plot) and M89, M59, M84, M60 and M87 at increasing radial velocities (bottom 
plot, M86 has v < 0). 



halo 



P2(R) 



oc 



3M 2 (fl) 
4irR 3 
M 2 (R)/M 2 
rf 



W [la^ — 77/(1 + 77)]} 



(7) 



where we assumed NFW density profile s and where the 
last equality is from lCole fc Lacevl l|l996|) . Since the mean 
density is a monotonously decreasing function of radius 



(i.e. the function of r\ in cq. \7\ decreases with 77), the 
most isolated halos will have the highest minimum values 
of 77. 

The output from GallCS includes the mass of each halo 
obtained from the FoF estimator, so that one can easily 
perform a first guess of which of the halos with mass close 
to that of Virgo are as isolated as appears the Virgo clus- 
ter. Table 3 lists for the 10 most massive halos in the 
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simulation, the parameters of its most perturbing neigh- 
bor (that with the lowest value of 77) using FoF quantities 
for radii and masses. Columns 2 to 4 list the virial radius, 
mass and circular velocity obtained using the more precise 
spherical overdensity (SO) method, which allow a direct 
comp arison with the value we obtained for the Virgo clus- 
ter in lMamon et alJ l|2004|) . 

Table 3. Most perturbing neighbors of 10 most massive 
halos in Gal ICS simulation 



add the Hubble flow. Therefore: 



rank 








R 




?7min 


(FoF) (Mpc) 


(1O 14 M ) 


(kms" 1 ) 


rf° 


M FoF 




1 


3.19 


16.7 


1503 


0.55 


0.011 


3.7 


2 


2.42 


7.4 


1143 


0.58 


0.003 


3.1 


3 


2.50 


8.1 


1180 


0.55 


0.020 


2.4 


4 


2.54 


8.5 


1200 


0.53 


0.011 


2.4 


5 


2.10 


4.8 


989 


0.62 


0.001 


2.9 


6 


2.00 


4.1 


942 


0.58 


0.001 


2.5 


7 


1.97 


3.9 


927 


0.57 


0.031 


3.4 


8 


1.78 


2.9 


842 


0.61 


0.002 


2.5 


9 


1.84 


3.2 


869 


0.58 


0.103 


2.2 


10 


2.03 


4.3 


958 


0.51 


0.003 


2.6 



As seen in Table 3, the most perturbing halo is always 
inside r\. Therefore, it is not clear that we can easily es- 
timate the exact virial radius and mass of each perturb- 
in g halo with th e same spherical overdensity method as 
in iMamon et alJ l|2004|) to produce a better estimate of 
isolated halos. As we expect rf° to be approximately pro- 
portional to rf oF , we take the results of Table 3 as good 
indicators of the degree of isolation of halos in the simu- 
lations. 

As shown in Table 3, the most isolated halo in the 
simulations is halo 1, the most massive halo. Picking the 
most massive halo also ensures that no other massive halos 
will distort the mock velocity-distance relation. 

3.3. Reseating the simulations to the scales of the 
Virgo cluster 

We now need to rescale the simulated clusters to Virgo. 
A reasonable way to put two halos in comparable units is 
to normalize distances with the virial radius, fiooj and 
velocities with the circular velocity of the halo at the 
virial radius, Uioo- We therefore convert from the sim- 
ulation frame to the Virgo cluster frame by rescaling 
the terms related to distance with the factor (omitting 
the '100' subscripts for clarity) ry/r s and the terms re- 
lated to velocities by the analogous factor vy/v s , where 
superscripts 's' and 'V refer to the simulation and the 
Virgo cluster, respectively. We adopt the Virgo cluster 
virial radius and circular velocity derived bv lMamon et all 
(2004) from X-ray observations ry = 1.65 h 2 , 3 Mpc, where 
h 2/3 = HoAee.Tkms^Mpc- 1 ), and vy = 780 kms" 1 . 

Moreover, the TV-body simulations give us comoving 
velocities u s of dark matter particles, to which we must 



D = D s — 

r s 

Vy 

V = U s 



H D 



(8) 
(9) 



where u s represents the comoving radial velocity of the 
particle in the simulation with respect to the observer in 
comoving units. 

Inverting equation ©, we need to place the observer 
at a distance 



Dobs = Dy — 
ry 



(10) 



to the halo used to mimic the Virgo cluster (in simula- 
tion units). Once we have D oos , we can select different 
observers by choosing different positions on a sphere of 
radius D oos centered on the halo and assigning to each 
position, the mean velocity of the 10 nearest halos to the 
observer, so as to give the observer an appropriate peculiar 
velocity. As our galaxies in the Virgo cluster extend up to 
fl max = 50 Mpc, we will need to use particles around the 
chosen dark matter halo up to a distance from the cluster 

Of Anax - Dy ~ 20ri 00 . 

It is interesting to check the peculiar velocities (or 
equivalently the systemic velocities) found for the simu- 
lated observer relative to the simulated cluster and com- 
pare to that of our Local Group of galaxies relative to 
the Virgo cluster, after rescaling the simulation radial and 
velocity separations to the Virgo cluster attributes as in 
equations JSJ and 10 • For 500 observers placed at a given 
distance from the halo, but in different random directions, 
we find a systemic velocity between the halo and the ob- 
server (i.e. the mean of 10 nearest halos) in Virgo units 
of 998 ± 118kms _1 (errors are the la dispersion) when 
we assume Dy = 16.8 Mpc (as suggested by the posi- 
tions of the early-type galaxies in Fig. ^ an d especially 
of M87, which lies at the center of the X-ray emission) 
and 1308 ± 155 kms -1 w hen we assume Dy — 21 Mpc, 
as in ISanchis et alJ ( 20021) . The values of the systemic ve- 
locity of the Virgo cluster in the literatu re range approx- 
imate ly between 900 and 1000 km s _1 ijTeerikorpi et alJ 
1992). Therefore, putting the Virgo cluster at 16.8 Mpc, as 
indicated by the early type galaxies, gives a better match 
to the systemic velocity of the cluster. In addition, if we 
calculate the number of all galaxies (early-type and late- 
type ones) per distance bin, we found the maximum of 
the distribution between 16 — 18 Mpc. Even if our sample 
lacks of completeness, this result, together with the pre- 
vious ones, lead us to choose Dy — 16.8 Mpc for the com- 
parison between the TV-body simulations and the Virgo 
cluster. 

3.4. Velocity-distance relation without distance errors 

We take all the objects lying within a cone of half-opening 
9° aligned with the axis connecting the observer to the 
center of the halo. 
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Because observed galaxy catalogs, such as the VCC 
ijBinggeli. Sandaee. k. TammanrJ [l985^ catalog of Virgo 
cluster galaxies, are magnitude-limited, we imposed upon 
our simulated catalog (of dark matter particles) the 
same magnitude lim it as the VCC, namely B < 18 
(|Binggeli et al.lll985|) . For this, we consider bins of dis- 
tance to the observer, and, in each distance bin, we re- 
move a fraction of dark matter particles equal to the pre- 
dicted fraction of particles fainter than B — 18 given the 
distance and the gaussia n-shaped luminosity function of 
spiral galaxies derived bv lSandage. Binggeli. k. TammanrJ 
l)l985|) . This method supposes that the luminosity func- 
tion is independent of environment (i.e. distance to the 
center of the Virgo cluster) and that those VCC galax- 
ies with redshift-independent distance measurements have 
the same magnitude limit as the VCC galaxies in general. 
This luminosity incompleteness is only noticeable at large 
distances, where it amounts to reducing the number of 
particles by less than 10%. 

Figure El shows the velocity-distance plot of the sim- 
ulated dark matter particles {small points), rescaled to 
Virgo units, and from which particles statistically fainter 
than the magnitude limit of the observed catalogs were re- 
moved. Superposed are the spiral galaxies {triangles) from 
ISolanes et alJ (|2002fl and the early- type galaxies {circles) 
compiled in Table El We have added to our s ample the 
spiral galaxies co mpiled bvlSolanes et alJ (2002) with reli- 
able Tully-Fisher (|Tullv & Fisherlll977l hereafter TF) dis- 
tances but no HI data. We list these galaxies in Table 4, 
with their coordinates, morphological type (from LEDA), 
systemic velo city and the final di stance with errors as- 
signed to it bv lSolanes et alJ ((2002) , who estimated errors 
in distance as the mean deviation between the rescaled 
mean distance and the different estimates used (with this 
method, galaxies with only one estimator cannot be as- 
signed an error and their distances are considered uncer- 
tain). We represent these objects as squares in Figure El 

Although FigureEJconcerns a single observer and a sin- 
gle halo to represent the Virgo cluster (i.e. halo 1, which is 
the most isolated), the results are very similar for all halos 
studied as well as for any observer whose velocity with re- 
spect to the center of the main halo is similar to the mean 
of the 500 random trials. In other words, FigureElprovides 
a realistic general velocity-distance relation for a Virgo- 
like cluster in a fla t Universe with a c osmological con- 
stant. Very recently. iKlvnin et al J <)2003h have produced a 
similar plot for dark matter particles from their adaptive 
mesh refinement simulation. Note th at the groups visible 
in the 3D phase space plot (Fig. 1 in lMamon et alJl2004[) 
are much less visible here and the projection places them 
at smaller distances from the center of the main cluster. 

Figure El indicates that the velocity field drawn by 
the observed spirals is noticeably different from the one 
followed by the dark matter particles of the simulation. 
Nevertheless, the figure shows that 1) the center of the 
halo corresponds to the location of greater galaxy density 
(not only for the early-type galaxies used to mark the cen- 
ter of the Virgo cluster), 2) the thickness of the halo in 



velocity space (the velocity dispersion) corresponds to the 
corresponding thickness of the central observed galaxies, 
and 3) the majority of the elliptical galaxies outside the 
cluster also appear to fit well the location of the particles 
in the outer infalling/expanding locus. These 3 facts sug- 
gest that both the distance to the halo and the rescaling 
factors are well chosen and that the differences between 
the Virgo cluster and the halo in the simulation should be 
attributable to errors arising from the distance estimates 
(as the velocity errors are negligible). It appears unlikely 
that the galaxies that, for a given distance, lie within more 
than 1000 km s -1 from the locus of dark matter particles 
belong to the low- or high-velocity tails of large velocity 
dispersion groups near the cluster, since in our simula- 
tions, there are no such high velocity dispersion groups 
outside the main cluster. 

3.5. Velocity-distance relation with distance errors 

Given that most of the galaxies outside the velocity field 
drawn by dark matter particles are spirals, we now build 
a new velocity-distance diagram of our simulated particles 
incorporating the mean relative distan ce error for the spi- 
ral galaxy sample. TSolanes et al.l l)2002j) combined different 
studies of the TF relation to calculate the distances of spi- 
rals in the Virgo cluster. These studies give dispersions of 
the TF relation of roughly 0.4 mag., corresponding to an 
uncertainty of 18% in relative distance. Figure |3| shows 
the corresponding velocity-distance plot for our simulated 
Virgo line of sight, with inclusion of gaussian relative dis- 
tance errors with oTnD = 0.20 (a slightly more conserva- 
tive value than 0.18). We omit particles beyond 40 Mpc, 
which corresponds to the distance (before folding in the 
distance errors) where the particles have a systemic veloc- 
ity approximately equal to the maximum of the observed 
galaxy sample. 

FigureEJindicates that the inclusion of distance errors, 
at the 20% relative rms level, allows us to reproduce fairly 
well the observed velocity-distance diagram. There are still 
several galaxies that lie far outside the general locus of the 
dark matter particles. Among them are two galaxies (VCC 
319, which is part of a galaxy pair and M90 = NGC 4569 
= Arp 76, which is tidally perturbed by a close neighbor) 
apparently in the foreground of the cluster between 7 and 
10 Mpc and with v < — 200 km s -1 , and four galaxies (IC 
3094, NGC 4299, which is part of a pair, NGC 4253 and 
VCC 1644) in the background between 34 and 42 Mpc 
and v < 800 km s -1 . Inspection of Figure El indicates that 
these four galaxies are at more than 2000 km s -1 from the 
expected locus of the Tolman-Bondi pattern given negligi- 
ble distance errors. It is highly unlikely that the peculiar 
velocities in interacting pairs can reach values as high as 
2000 kms -1 , hence the distances to these galaxies appear 
to be grossly incorrect. 

We list the most discrepant galaxies in Table 5. Most 
of these galaxies have a single distance measurement (the 
elliptical VCC 319 has a single Sersic fit distance esti- 
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Table 4. Spiral galaxies without Hi data 



Galaxy 


RA 


Dec 


T 


V 


D 


VCC 


NGC 


(J2000) 




(km s" 1 ) 


(Mpc) 


132 




12 h 15 m 03 ! :9 


+13°01'55" 


8 


1965 


10.76t£S 


343 




12 h 19 m 22 ! :0 


+07° 52' 16" 


8 


2336 


23.77 


912 


4413 


12 h 26 m 32 ! :2 


+12°36'39" 


2 


-13 


15. 63^2. 20 


952 




12 h 26 m 5S?5 


+09°52'58" 


5 


853 


4-77i«° 


1605 




12 h 35 m 14 ! ;5 


+10°25'53" 


5 


951 


24.66t 6 5 ;^ 


1673 


4567 


12 h 36 m 32=9 


+11°15'28" 


4 


2133 


26.30 


1933 




12 h 42 m 44 s 7 


+07°20'16" 


2 


2271 


35.16t ; 32 
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Fig. 2. Simulated and observed velocity-distance diagrams. The observer is at (0, 0) and sees a cone of angular radius 
9°. Dots represent the velocity field traced by the particles (1 in 50 for clarity) in the cosmological A^-body simulation 
in the direction of the most massive halo, for which distances and velocities have been rescaled (eqs. [H] and [5]) 
incorporating a Hubble flow with Ho = 66.7kms _1 Mpc" 1 , as in the simulation. Superposed are galaxies of the Virgo 
cluster. Circles and triangles represent early-type and late-type galaxies, respectively. The size of the triangles is 
proportional to the H i-deficiency of the spiral galaxies measured in units of the mean standard deviation for field 
objects (= 0.24). Squares represent additional spiral galaxies with no H i-deficiency data. Open symbols indicate 
galaxies with uncertain distances (with distance errors greater than 5 Mpc or spirals with only one distance estimate) . 
Dashed lines show the unperturbed Hubble flow with Hq = 66.7 and 70kms _1 Mpc -1 , respectively (going upwards). 
The locations of the vertical spikes of simulated particles away from the Virgo cluster (assumed to lie at 16.8 Mpc) 
are not meaningful and should not be compared to the measured galaxy distances. 



mator, which is known to be inaccurate — see Table 1). 
However, it is a surprise to find a bright object such as 
NGC 4569 (M90) in this list, especially since it has 6 differ- 
ent Tully-Fisher distance measurements, 5 o f which yield a 
dista nce between 7.9 and 10.2 Mpc and one ijGavazzi et alJ 
119991) giving D = 16.1 Mpc, which is what we expect from 
Figure |2 

Since the 20% relative errors correctly reproduce the 
envelope of particles from the cosmological simulation, in- 
creasing the relative rms error to say 30% will not explain 
the distances to these discrepant spirals. We therefore con- 



clude that the distances of the galaxies listed in Table 5 
must be highly inaccurate. 

In their dynamical, Tolman-B ondi. modeling of the 



Virgo region, I Sanchis et alJ ( 2002h only used those spiral 
galaxies believed to trace the outer envelope of galaxies 
expanding out of the Virgo clust er. Figure El clearly in- 
dicates that the galaxies used bv ISanchis et alJ (roughly 
aligned between distances of 34 and 45 Mpc) lie at 
typically 600 km s -1 lower velocities for their distance 
than predicted by our simulated velocity-distance model. 
Figure [3 suggests that the alignment of these galaxies in 
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Fig. 3. Same as Figure 2, but incorporating gaussian relative distance errors of ah\D = 0.2 for the dark matter 
particles. 



Table 5. Grossly incorrect galaxy distances 



Galaxy 


RA Dec 


T 


V 


D 


VCC NGC 


(J2000) 




(kms- 1 ) 


Mpc 


213 


12 h 16 m 56!0 +13° 37' 33" 


5 


-278 


40.9 


319 


12 h 19 m 0r!7 +13°58'56" 


-4 


-256 


7.6 


491 4299 


12 h 21 m 40 s 9 +11 ° 30 '03" 


8 


112 


34.8 


1524 4523 


12 h 33 m 48!0 +15°10'05" 


9 


158 


36.1 


1644 


12 h 35 m 51!8 +13°51'33" 


9 


646 


41.5 


1690 4569 


12 h 36 m 49:8 +13°09'48" 


2 


-328 


9.4 



VCC 1690 = NGC 4569 = M90 = Arp 76. 



the velocity-distance plot may be fortuitous, given the 
~ 20% relative dist ance errors. T his fortuitous alignment 
forced the model of lSanchis et al.l to shift the center of the 
Virgo cluster to greater distances. 

3.6. Tolman-Bondi estimates of the Virgo cluster mass 
and distance 

We now use the cosmological iV-body simulations to quan- 
tify the errors on the distance and mass to the Virgo 
cluster obtained by fitting a Tolman-Bondi model to the 
velocity-distance data outside the cluster core. We con- 
struct simulated pseudo-galaxy samples by randomly se- 
lecting as many points from the cosmological iV-body sim- 
ulation as there are galaxies in the Virgo cluster sample. 
We apply fixed relative distance errors in the range 0%- 
30% and, for each relative distance error, we construct 
50 realizations of the Virgo sample. We identify pseudo- 
galaxies presumably belonging to the core by building his- 



tograms of the number of pseudo-galaxies per 2 Mpc dis- 
tance bin. Such histograms peak at the core of the halo, 
in a region of width ~ 4 Mpc for the no error set and up 
to ~ 12 Mpc for the 30% error set. We avoid the virialized 
core by only considering pseudo-galaxies lying on the back 
side of the identified cluster core. 

In the Tolman-Bondi model used by ISanchis et all 
galaxies are treated as test particles moving in a point- 
mass gravitational potential that is independent of time. 
The model velocity-distance curve is a function of the age 
of the Universe, to, the distance to the cluster, Dy, the 
point mass of the cluster, My, and the systemic velocity 
of the cluster relative to the Local Group of galaxies, vy. 
Only three of these four parameters are independent, but 
given the poor current velocity-distance data, it is bet- 
ter to use only two free parameters and fix the third one. 
Following lSanchis et all we have fixed vy to the mean ve- 
locity of the main halo for the 50 realizations (998 km s _1 
for a distance of 16.8 Mpc, see section 3) and we have 
allowed My and Dy to vary freely. The resulting to will 
help us to see the viability of the results. 

Table 6 shows the results of the fits to the po int-mass 
Tolman-Bondi relation (see I Sanchis et all [2002) for the 
different sets. For zero distance errors, the model repro- 
duces fairly well t he input data of D y = 16.8 Mpc and 
My = 10 14 - 35 M Q |Mamon et, alJl2004 - as well as the age 
of the Universe: to = 13.7 Gyr from the GallCS cosmologi- 
cal simulations (which happens to perfec tly match the age 
recently obtained bv lSpergel et alJl20o3 from the WMAP 
CMB experiment). The Table also clearly indicates that 
the errors in distance contribute to a slightly larger dis- 
tance to the cluster (by up to 2%), combined with a much 
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Table 6. Best fitting parameters for each dataset 



Error 


D v 




log(Mv) 


to 




(Mpc) 




1 
I 


(Gyr) 




mean 


a 


mean 


a 




0% 


16.7 


0.5 


14.5 


0.2 


14.1 


5% 


16.8 


0.5 


14.6 


0.3 


13.9 


10% 


16.9 


0.6 


14.9 


0.4 


12.8 


15% 


17.0 


0.8 


15.1 


0.5 


12.1 


20% 


17.1 


0.9 


15.5 


0.5 


10.1 


25% 


17.1 


1.0 


15.7 


0.4 


9.0 


30% 


17.1 


1.0 


16.1 


0.5 


6.7 



larger mass for the cluster (up to a factor 40!). This dra- 
matic rise in the mass decreases the resulting to to unac- 
ceptably low values. Therefore, distance errors can cause 
the fitting of Tolman-Bondi solutions to the velocity vs. 
distance of galaxies to yield incorrect cluster parameters. 

For the case of real data on the Virgo cluster, the 
same fitting procedure leads to My ~ 3 x 10 16 M Q and 
Dy ~ 17.4 Mpc, which imply £ = 5.2 Gyr. Note that 
the present fitting procedure takes into account all galax- 
ies beyond th e identified core of the cluster, whereas 
ISanchis et alJ only used those galaxies that are aligned be- 
tween distances of 34 and 45 Mpc in the velocity-distance 
plot, deducing a very different distance to the Virgo clus- 
ter (between 20 and 21 Mpc). The comparison of the val- 
ues of to from the real data and from the simulated data 
convolved with different relative distance errors, suggests 
that the TF distances may have up to 30% relative errors. 

As shown above, the estimation of the distance to the 
cluster provided by the point mass model is relatively good 
even for appreciable errors in relative distances. Therefore, 
one could still obtain a reasonable value for the mass of the 
cluster by taking Dy fro m the fit (i.e., 17.4 Mpc ), together 
with vy = 980 km s" 1 iTeerikorpi et alJll99^) and to = 
13.7 Gyr ijSpergel et al.ll2003|) as input parameters for the 
model and computing My as a result. In this way, we 
obtain a value of My = 7 x 10 14 M Q . Notice that the 
effective mass of the point mass model does not have to 
coincide with the virial mass of the cluster, as the latter 
one represents only the virialized part of the cluster, while 
the former gives us an idea of the dynamical mass of the 
whole Virgo region. 

4. Discussion 

As mentioned in Sec. 1 and clearly seen in Figure[21 there 
is a non negligible number of H I-deficient galaxies distant 
from the cluster core, in particular at distances of 10 and 
28 Mpc from the Local Group. We now examine different 
explanations for the presence of these H I-deficient galaxies 
outside the core of Virgo: 

1. very inaccurate distances (perhaps from the inade- 
quacy of the Tully Fisher relation for galaxies under- 
going ram pressure or tidal stripping) for galaxies that 



are in fact within the Virgo cluster (in the core or 
more probably in a rebound orbit) and which would 
have lost some of their gas by ram pressure stripping 
from the intracluster gas; 

2. tidal interactions from members of a same group; 

3. tidal interactions from neighbors; 

4. recent minor mergers heating the gas; 

5. overestimated H I-deficiency in early- type spirals. 

In Table 7 are listed a ll galax ies from the 
compilation of ISolanes et al ] ll2002l) which sat- 
isfy at least one of the following criteria: 

• D < 13 Mpc , 

• 21< D < 50 Mpc , 

• o/e 10 o > 1 , 

where 9 is the angular distance to M87 and 
0ioo = r ioo/-Dv = 5? 6 is the angular virial radius 
of the Virgo cluster. Columns (1) and (2) give the galaxy 
names, the coordinates in columns (3) and (4), the galaxy 
type, total blue extinction corrected magnitude and 
optical diameter (all three from LEDA) are respectively 
in columns (5), (6), and (7), the projected distance to 
M87 in units of the virial radius (rioo) of the cluster in 
column (8) , the velocity relative to the Local Group taken 
from the Arecibo General Catalog (a private dataset 
maintained by R. Giovanelli and M. P. Haynes at the 
University of Cornell), in column (9) and co lumn (10) 
gives the adopted distance ^Solanes et alJl2002^ . The last 
3 columns of Table 7 are explained in Sect. 14.11 

4.1. Are the outlying H I-deficient galaxies truly too far 
from the Virgo cluster to have had their gas ram 
pressure stripped? 

If one wishes to explain the outlying H I-deficient galaxies 
as galaxies that have lost their interstellar gas by its strip- 
ping by ram pressure from the hot intracluster gas, then 
these galaxies must have passed through the core of the 
Virgo cluster where the intracluster gas is dense and the 
gal axy velocities l arge so that the ram pressure is largest. 

iMamon et alJ l|2004f) found that objects that have 
passed through the core of a structure in the past can- 
not be at distances greater than 1-2.5 virial radii from 
the cluster today. With their estimate of the distance (16.8 
Mpc) and the virial radius (1.65 Mpc) of Virgo, this means 
that galaxies that have passed through the core of Virgo 
cannot lie at a distance from the Local Group greater than 
18.5 or 20.9 Mpc (for r re b/^ioo = 1 or 2.5, respectively) 
nor can they lie closer than 15.1 or 12.7 Mpc from the 
Local Group. 

For the 3 H I-deficient spirals at 28 Mpc to have passed 
through the core of Virgo, one would require that their dis- 
tances be each overestimated by 52% (34%), which corre- 
sponds to 2.3 a (1.6 a) events for their distance moduli, for 
^rebAioo = 1 (2-5), respectively. Such errors are possible, 
as shown in Sect. 13.51 Note that none of the distances of 
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Galaxy 


RA 


Dec 


T 




£>25 


e 


V 


D 


N 


Pi 


P2.5 


VCC 


NGC 


(J2000) 










(kms" 1 ) 


(Mpc) 










4064 


12 h 04 m ll!8 


18°26'33" 


1 


11.73 


4f0 


8? 8 


837 


9.8 


942 


0.00 


0.81 


522 


4305 


12 h 22 m 03!5 


12°44'27" 


1 


12.87 


2'0 


2? 2 


1814 


40.7 


200 


0.00 


0.00 


524 


4307 


12 h 22 m 06!3 


09°02'27" 


3 


11.84 


3.'5 


4?0 


913 


27.4 


1846 


0.41 


0.48 


559 


4312 


12 h 22 m 32!0 


15°32'20" 


2 


11.76 


4.'7 


3? 7 


47 


10.8 


754 


0.95 


1.00 


713 


4356 


12 h 24 m 15 ! ;9 


08°32'10" 


6 


13.02 


2'6 


4? 2 


998 


29.4 


1377 


0.19 


0.23 


979 


4424 


12h27 m 13!3 


09°25'13" 


1 


11.99 


3.'4 


3?1 


314 


4.1 


106 


0.00 


0.00 


1043 


4438 


12 h 27 m 46 s 3 


13°00'30" 


1 


10.55 


8.' 7 


1?0 


-45 


10.4 


755 


1.00 


1.00 


1330 


4492 


12 h 30 m 5g!9 


08°04'41" 


1 


13.04 


1.'9 
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the 13 galaxies listed in Table 7 are based upon Cepheid 
measurements, which are much more precise than Tully- 
Fisher estimates. On the other hand, given the presence in 
our sample of some spirals with grossly incorrect distances 
(see Sect. one may wonder if the distances to the 3 
spirals at 28 Mpc could also be grossly incorrect: they 
could lie at the cluster distance of 16.8 Mpc, thus leading 
in this case to distances overestimated by 66%, which for 
20% rms relative errors correspond to three 2.8 a events for 
their distance moduli, which appears unlikely, unless the 
Tully-Fisher error distribution has non-gaussian wings. 

Similarly, the three foreground H I-deficient spirals at 
~ 10 Mpc could have bounced out of the core of the Virgo 
cluster, as they lie on the outer edge of the envelope of 
the particles after inclusion of 20% rms distance errors in 
Figurc|3| Indeed, were they bouncing out of the Virgo clus- 
ter in the foreground at 15.1 (12.7) Mpc, one would then 
have three cases of 34% (21%) errors, each corresponding 
to —2.3cr (— 1.3<7 ) events for their distance moduli, for 
frebAioo = 1 (2-5), respectively. Alternatively, they could 
lie in the cluster proper, at 16.8 Mpc, which would cor- 
respond to three — 2.8 <r events for their distance moduli, 
which again appears unlikely, unless the Tully-Fisher error 
distribution has non-gaussian wings. 

The offset of a galaxy relative to the Tolman-Bondi 
locus in Figure can be caused by a large distance error 
or a large peculiar velocity relative to the general flow 
at that distance, or a combination of both. To see which 
effect is more important, one can consider the N simulated 
particles with distance (with errors folded in) and velocity 
respectively within 3 Mpc and 200 km s -1 of each galaxy 
in Table 7, as well as angular distance to the halo center 
within 1?5 of that of the galaxy relative to M87, and asked 
what fractions Pi and P2.5 of these particles were located 
(before the distance errors were folded in) within 1 or 2.5 
times the virial radius of the halo, respectively. 

These fractions must be taken with caution, because 
the groups falling in or moving out of the Virgo cluster are 



at different distances and with different relative masses be- 
tween the simulations and the observations. For example, 
a true group of galaxies in the background of the Virgo 
cluster will end up with too high values of Pi and P2.5 if 
there is no group at the same distance in the simulations, 
while a Virgo cluster galaxy with a measured distance of 
10 Mpc or so beyond the cluster will have too low val- 
ues of Pi and P2.5 if the simulations include a group at a 
distance close to this measured distance. Nevertheless, Pi 
and P2.5 should provide interesting first order constraints. 

Of course, this method of locating galaxies, implies 
that galaxies NGC4064, NGC 4580 and NGC 4586, at 
angular distances from M87 beyond 9±oo + 1?5 = 7?1 have 
Pi = 0. Interestingly, all three galaxies are likely to be 
within 2.5rioo (i.e. P2 is close to unity). 

Among the 5 foreground galaxies with 9 < 5?6, NGC 
4312, NGC 4438, M90 and NGC 4606 are all likely to be in 
the Virgo cluster and all very likely to be close enough that 
they may have passed through its core. Only NGC 4424, 
apparently at distance 4.1 Mpc, is too close to actually be 
located within 2.5 virial radii from the Virgo cluster. 

Among the 5 background galaxies with 9 < 5? 6, all 
but NGC 4305, lying at 40 Mpc, have a non negligible 
probability of having crossed the cluster. Only VCC 1569 
measured at 23 Mpc is highly likely to be located within 
the virial radius of the Virgo cluster. Interestingly, the 3 
galaxies at measured distances around 28 Mpc have very 
different values of Pi and P2.5. The closest one, NGC 4307, 
has one chance in two of being located near the cluster and 
40% probability of being within the cluster itself. However, 
the furthest one, NGC 4356, has one chance in 5 of belong- 
ing to the Virgo cluster and less than one chance in four 
of being within 2.5 virial radii. The third one, NGC 4492, 
has only one chance in 9 of being within 2.5 virial radii of 
the Virgo cluster. This low probability is not surprising, 
given that it is located close to the infall/expansion zone 
of Figure □ 
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4.2. Are the outlying H l-deficient galaxies located 
within small groups? 

The alternative to ram pressure stripping of gas is removal 
of gas by tidal effects. One can envision various scenarios: 
1) tidal stripping of gas beyond the optical radius; 2) tidal 
compression of gas clouds leading to gas transformation 
into stars; 3) tidal stripping of the gas reservoir infalling 
into the disk. Such tides can occur from interacting galax- 
ies or from the cluster or group itself. We begin by asking 
whether the Hi-deficient galaxies apparently outside of 
Virgo lie in groups. 

Figure 01 shows the environment of the H I-deficient 
galaxies in three equal logarithm distance bins, ±1.28 x 
crlnD wide, centered around 10 Mpc (the apparent dis- 
tance to the foreground H I-deficient galaxies), 16.8 Mpc 
(the distance to the Virgo cluster) and 28 Mpc (the ap- 
parent distance to the background H I-deficient galax- 
ies), as well as in two radial velocity bins. The log dis- 
tance (velocity) bins are chosen wide enough that if the 
parent distribution function of the log distance (veloc- 
ity) errors were gaussian, one would have a probability 
2 erf (1.28) — 1 = 0.86 of having a galaxy with a true dis- 
tance (velocity) within the interval appearing to lie at the 
center of the interval. 

Among the 5 foreground Hl-dcficient spirals, 4 (NGC 
4312, NGC 4438, and M90, which, as noted in Sect. ESI 
is far off the Tolman-Bondi relation, and NGC 4606), ap- 
pear superposed (within dotted circle in upper-left plot) 
with the galaxies at the distance of the main structure of 
the Virgo cluster, which suggests that their distances may 
be seriously underestimated, unless their happens to be a 
group of galaxies precisely along the line of sight to the 
center of Virgo lying at 10 Mpc. 

One may be tempted to identify the 3 H I-deficient spi- 
rals close to M87 (NGC 4438, M90, and NGC 4606) with a 
group of ellipticals seen in the upper- left plot of Figure 0] 
The velocity dispersion of this group would explain why 
these 3 galaxies are so far off the Tolman-Bondi solution 
given in Figure However, none of the ellipticals in this 
region have accurate (< 38%) relative distance estimates 
(cf. Table 2) and some or all may actually belong to the 
Virgo cluster. 

Interestingly, 3 of the 4 H I-deficient spirals at 28 Mpc 
(NGC 4307, NGC 4356 and NGC 4492) lie in a dense re- 
gion, extending to the SSW of the position of M87 (lower- 
left plot), which is also seen in the middle distance bin 
(middle-left plot) corresponding to the distance to M87, 
but with 3 times fewer galaxies. Although the eastern edge 
of the region, and particularly the H I-deficient NGC 4492, 
coincide with the position of M49 (NGC 4472), it appears 
distinct from M49, which roughly lies at the middle of a 
group of galaxies (middle-left plot) that does not seem to 
be associated with the dense region at 28 Mpc (lower-left 
plot). 

We now compare the Tully-Fisher distance estimator 
with the rcdshift distance estimator D v = v/Hq. If a struc- 
ture outside of the Virgo cluster lies at a distance D within 



a region with ID velocity dispersion <j v , the error on the 
distance estimator D v = v/Hq will be set by the disper- 
sion of the peculiar velocities: 

SD V = <j v /H . (11) 

The redshift distance estimator D v will be more accurate 
than the direct distance estimator when 

HoD> sSjD ' < 12 > 

where 8D /D is the relative error of the direct distance esti- 
mator. For Tully-Fisher distance measurements, SD/D ~ 
20% and a field velocity dispersi on o f 40 — 85kms -1 (as 
obser ved by lEkholm et al.l l200l| and iKarachentsev et alJ 
120031 as well as predicted bv lKlvpin et al.l 12003^1 . eaua- 
tion (|12|) yields a maximum distance for believing individ- 
ual direct distances of 3 to 7 h^,^ Mpc, so that individual 
galaxy distances of 28 Mpc are much less reliable. 

This change in optimal distance estimators is illus- 
trated in Figure n where only one of the background H I- 
dcficient galaxies has a radial velocity much larger than 
the mean of the Virgo cluster. The two plots on the right of 
Figure 0] show that the structure to the South- Southwest 
(SSW) of the Virgo cluster lies close to Virgo in velocity 
space (i.e. the concentration is in the right middle plot). 
In redshift space, the 3 background H I-deficient galaxies 
to the SSW of M87 remain in a similar concentration of 
galaxies (middle right plot of Fig. 0} as in distance space 
(lower left plot), but this concentration is at similar radial 
velocities as the Virgo cluster. However, only one (NGC 
4307) has a reasonably high chance (Pi in Table 7) of be- 
ing located in or near the Virgo cluster, while one (NGC 
4356) has less than one chance in 4 of being in or near 
the cluster (from the analysis of Sect. 14. l|) . and one (NGC 
4492) has a radial velocity indicative of a background ob- 
ject expanding away from the Virgo cluster (this galaxy 
shows up as the empty triangle in Fig. |3 very close to the 
Tolman-Bondi locus) , and has only roughly one chance in 
9 of being even close to the cluster. But the analysis of 
Pi and P2.5 °f Sect. 14. ll should be superior to the analy- 
sis of Fig. 01 because it takes into account the correlation 
of distance and velocity. In summary, the 3 background 
H I-deficient galaxies to the SSW of M87 may be suffering 
from tidal effects of a background group of galaxies. 

4.3. Close companions and mergers 

According to a search in SIM BAD, we find that (at least) 4 
of the 13 H I-deficient galaxies in Table 7 have companions 
whose mass and proximity may generate tides that could 
remove some of the neutral Hydrogen, either by stripping 
the gas beyond the optical radius (without stripping the 
stars, thus leaving an H I-deficiency) or by compressing 
the diffuse neutral hydrogen, converting atomic gas into 
molecular gas, then into stars. Of course, the presence of 
a companion in projection, even with a similar velocity, 
does not guarantee that this companion is in physical in- 
teraction with the galaxy in question, as projection effects 
are important in clusters. 
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Fig. 4. Positions of galaxies within 9° of M87 in 3 distance (left) and velocity (right) bins. Ellipticals and spirals are 
shown as circles and triangles, respectively. Messier ellipticals are shown as filled circles, while M49 (NGC 4472) is 
also highlighted as a thick dotted circle (bottom plots, in green in the Electronic version of the journal). Hi-deficient 
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Table 8. Companions to H i-deficicnt galaxies in the out- 
skirts of Virgo 



Galaxy 


Companion 


Aw 


Am c T 


d/r opt 


M 


K 


VCC NGC 


VCC NGC 


(kms" 1 ) 










522 4305 


523 4306 


-380 


O.T 


2.8 


R? 




1043 4438 b 


1030 4435 


-214 


1.0 


1.0 


I 


I 


1690 4569 c 


1686 d 


1354 


4.2 


1.1 


R 


R 


1859 4606 


1868 4607 


593 


0.6 


2.7 


I? 





Notes: a B-band; 6 also Arp 120; c also M90 and Arp 76; d also 
IC 3583. 



Table 8 shows the properties of the companions to 
these 4 galaxies. Column 5 is the difference in velocity 
relative to the Local Group; column 6 is the difference 
in extinction-corrected total /-band magnitude; column 7 
is the angular distance between the two components, in 
units of the optical radius (D25/2) of the major member; 
columns 8 and 9 are indicat ors for morphological and kine- 
matical llRubin et alJll99^) disturbances (R — > regular, I 
— * irregular) for the major component. Table 8 indicates 
that NGC 4438 has a very close major companion, which 
is likely to be responsible for its irregular morphology and 
internal kinematics. NGC 4305 and NGC 4606 have mod- 
erately close major companions. M90 has a close very mi- 
nor (m 1/50 in mass) companion. 

Moreover, NGC 4307 and NGC 4356 are close in 4D 
space: their projected separation (at the distance of M87) 
is 215 kpc and their velocity difference is 85 kms , so 
they may constitute an interacting pair. Similarly, NGC 
4580 and NGC 4586 are at a projected separation of 311 
kpc with a velocity difference of 254 kms -1 , so they too 
may constitute an interacting pair. 

Finally, two othe r galaxies in Table 7 have unusually 
low inner r otation fe'ubin et al lll999T) : NGC 4064 and 
NGC 4424. iKennev et all l|l996l) argue that NGC 4424 
has undergone a merger with a galaxy 2 to 10 times less 
massive and presumably the same can be said for NGC 
4064. 

We can obtain some clues about the origin of the H I 
deficiency of the outlying galaxies by studying their star 
formation rates. We adopt the, obviously simplified, point 
of view that ti dal interactions are beli eved to enhance star 
formation (e.g. lLiu fc Kennicuttll995|) . while ram pressure 
stripping should reduce star formation. We therefore verify 
the tidal interaction hypothesis by checking for enhanced 
star formation rates. Because tidal interactions operate 
on scales of order the orbital times of interacting pairs, 
typically 250 Myr or more, we do not use the Ha star 
formation rate indicator, as it is only sensitive to very 
recent star formation, but focus instead on the broad-band 
color B—H. We have taken B and H magn itudes from the 
GOLDMine database (|Gavazzi et alJl2003|) . 

Figurc[5]shows the B—H colors as a function of if-band 
luminosity for the 12 NGC galaxies of Table 7 (we have 
omitted VCC 1569, which is too undcrluminous to be in 
GOLDMine) and for other Virgo galaxies. The colors and 



luminosities in this figure were corrected for absorption 
using internal and galactic ijSchlegel. Finkbeiner. fe Davisl 
1998) absorption obtained fr om LEDA, assuming selective 
extinction coefficients from ICardelli. Clavton. fc M athis 
( 1989). Figure[3clearly displays a general color-luminosity 
relation, i.e. the trend that more luminous galaxies are 
redder. Comparing galaxies of similar mophological types, 
one checks that among the four candidates for possess- 
ing tidal companions, NGC 4305, NGC 4438, NGC 4569 
(M90), and NGC 4606 (see Table 8), the first three are 
blue for their luminosities, assuming their Tully-Fishcr 
distances, but given that the latter three might well lie 
in Virgo (see column Pi in Table 7), we also check that 
at the distance to the Virgo cluster these 3 galaxies are 
very blue (NGC 4438), extremely blue (M90) and just 
blue (NGC 4606). The blue colors for these 4 galaxies, 
given their luminosity and morphological type, confirms 
the tidal hypothesis for all four objects. 

Interestingly, the possibly interacting galaxies 
NGC 4307 and NGC 4356 have blue colors confirming 
their interaction if they are indeed in the background 
of Virgo as their Tully-Fisher distances suggest, but 
have red colors as expected for ram pressure stripped 
galaxies if they lie in the Virgo cluster. The values of 
Pi in Table 7 cannot distinguish between these two 
possibilities. The situation is even more uncertain for the 
possibly interacting galaxies NGC 4580 and NGC 4586, 
whose Tully-Fisher distances place them slightly behind 
and in front of the Virgo cluster, respectively. If they both 
lie at the distance of the Virgo cluster, the first one would 
have a normal color, while NGC 4586 would be blue as 
expected if it is indeed interacting with NGC 4580. 

The situation is also confusing for the two galaxies 
with possible recent minor merging: NGC 4064 appears 
normal or red if it lies in the foreground as suggested by 
its Tully-Fisher distance, but blue if it lies at the distance 
of the Virgo cluster, and NGC 4424 appears normal or red 
if it lies as close as its Tully-Fisher distance suggests but 
quite blue if it lies at the distance of Virgo. One would 
expect the minor merger to generate a blue nucleus if the 
merging took place recently enough, but this may have 
little effect on the overall color of the galaxy. 

Moreover, one may contest that tides will remove stars 
as well as gas from the disk of a spiral galaxy. However, the 
estimates of H I-deficiency presented here are based upon 
single-dish (Arecibo) radio observations, which often do 
not have sufficient angular resolution to map the gas (the 
half-power beam of Arecibo subtends 3'3, close to the typ- 
ical optical angular diameters of the H I-deficient galaxies 
- see Table 7). One can therefore easily imagine that 
fairly weak tidal effects will deplete the neutral Hydrogen 
gas beyond the optical radius without affecting the stars 
in the spiral disk, thus leading to an H I-deficiency when 
normalized relative to the size of the optical disk. 

Two galaxies among our H I-dcficicnt Virgo outliers 
list ed in Table 7 (NGC 4 438 and M90) have been studied 
bv ICavatte et all l|l994l) . who computed H I-deficiencies 
of Virgo spirals, based upon high resolution 21cm VLA 
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Fig. 5. Color-luminosity diagram for the 12 of the 13 H i-deficient galaxies of Table 7 (VCC 1569 was excluded because 
it is not bright enough to be included in GOLDMine), possibly outside the Virgo cluster (large filled symbols) or placed 
at the distance of the Virgo cluster (large open symbols). Other Virgo galaxies are shown as small open symbols. 
Triangles, squares and pentagons (respectively in red, green, and blue in the Electronic version of the jo urnal) refer to 
galax ies of morphological type 1 — 2, 3 — 4 and 5 — 6, respectively. The data is taken from the GOLDMine ijGavazzi et alJ 
2003) database, and dereddend e d usin g the internal and Galactic extinctions from LEDA with the selective extinctions 
of lCardelli. Clavton. fc Mathisl l|l989h . 



observations. Both galaxies present very small 21cm to 
optical diameter ratios, indicating that the gas deficiency 
sets in at small radii and cannot be caused by tidal effects. 
Additional high- resolution 21 cm observations at the VLA 
of the other galaxies in Table 7 will obviously allow con- 
firmation of this result. 



4.4. Are the outlying H l-deficient galaxies truly 
gas-deficient? 

The H I-defic iency estimator has been defined by several 
authors fe.g.lCha maraux et al.ll980tlHavnes fc Giovanellil 
ll984tlSolanes et al.ll996|) as the hydrogen mass relative to 
the mean for galaxies of the same morphological type and 
optical diameter: 



DEFi = (log M m (D opt ,T)) - logM m , 



(13) 
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where Mht is the H I mass of the galaxy in solar units. 
ISolanes et al.l l)2002j) applied a different H I-deficiency es- 
timator based on the hybrid mean surface brightness ex- 
pected for each galaxy type independently of its diameter, 
which has the advantage of being independent of the dis- 
tance to the galaxies: 



DEF 2 = (logEHi(T))-logE H i , 



(14) 



where Shi is the hybrid mean surface brightness, defined 
as 



Shi = 



7 opt 



(15) 



where -Fhi and 9 opt are the HI flux and optical angular 
diameter, respectively. 

Could the seemingly Hl-deficient galaxies in the fore- 
ground and background of Virgo have normal gas content 
according to equation l | 1 3l) but appear deficient with equa- 
tion (|14|> ? ISolanes et alJ checked the statistical agreement 
of these two H I-deficiency estimators by fitting a linear 
model between the two parameters and found a regres- 
sion line with slope and intercept close to the unit and 
zero respectively, with a scatter or 0.126. Figure El shows 
this correlation, where those galaxies possible located out- 
side the Virgo cluster are shown as solid circles. The plot 
clearly indicates that the outliers are also deficient with 
the DEFi definition of H I-deficiency. 




Fig. 6. Comparison of H i-deficiency estimators. Galaxies 
possibly located outside the Virgo cluster are shown as 
filled circles. The vertical and horizontal dashed lines 
shows the 3er cutoff for DEFi and DEF2, respectively. 

However, Table 7 clearly shows that 7 out of our 
13 Hl-deficient possible outliers appear to be Sa galax- 
ies (T = 1), some of which might be misclassified 



SO/a or even SO galaxies, which are known to have 
little cold gas. The stronger H I-deficiency of early- 
type spi rals in the Virgo cluster has been alrea dy no- 
ticed by iGuiderdoni fc Rocca-Volmerangel l)l985h . Some 
of these Sa galaxies may not be deficient relative to 
field SO/a or SO galaxies. Indeed, the mo st recent study 
l)Bettoni. Galletta. fc Garci'a-Burilloll2003j) of the HI con- 
tent of galaxies as a function of morphological type in- 
dicates that the logarithm of the HI mass normalized to 
luminosity or to square optical diameter decreases increas- 
ingly faster for earlier galaxy types, therefore its slope for 
Sa's is considerably larger than for Sc's, and hence SO/a 
galaxies misclassified as Sa's will appear more H I-deficient 
than later-type galaxies misclassified by one-half of a mor- 
phological type within the Hubble sequence. 

Moreover, independent estimates of morphological 
types by experts led to a (T) = 1.5 (rms) fo r galaxies with 
optical diameters > 2' l|Naim et al.lll995|) . The galaxies 
listed in Table 7 have optical diameters ranging from 2' 
to 10' plus one (VCC 1569) with D = 0'.8. There are 
two galaxies (NGC 4305 and NGC 4492) with T = 1 and 
D < 2', i.e. where the uncertainty on T is probably as 
high as 1.5. 

For example, NGC 4424 and NGC 4492 are marginally 
H I-dcficicnt and are prime candidates for losing their sta- 
tus of gas-deficient if they are misclassified lenticulars (by 
AT = 1, the typical scatter given by LEDA). 

We must thus resort to visualizing the possibly mis- 
classified galaxies. Figure [7| displays the snapshots of the 
7 Hl-deficient outliers classified as Sa (T = 1). We now 
discuss the morphologies of these 7 galaxies. 

NGC 4064 This galaxy has a very large high surface 
brightness inner disk or flattened bulge with little trace 
of spiral arms. NGC 4064 may lie between SO and SO/a. 

NGC 4305 This galaxy is of type Sa or later. It has a 
nearby major companion. 

NGC 4424 This galaxy has an importan t bar and the 
size o f the bulge is difficult to estimate l|Kennev et alJ 
1996 fit the bulge-disk ratio to bet ween 0.2 and .6). It 
shows little traces of spiral arms. iKennev et alJ write 
"NGC 4424 is classified as an Sa rather than an SO 
because the outer stellar disk is not structureless". It 
may well be an SO/a. 

NGC 4438 This galaxy, with a large bulge, is tidally 
distorted by a com panion just North of the frame. 
iGavazzi et al.1 1 2000h fit a pure r 1 / 4 law to this galaxy, 
suggesting that its morphology is of an earlier type 
than Sa. 

NGC 4492 This galaxy is viewed face-on, hence the disk 
is not very visible. The bulge is relatively large and 
their are signs of a weak spi ral arm. NGC 4492 m ay 
lie between SO/a and Sa, but IGavazzi et al.1 l|200fjft fit 
a bulge to disk ratio of 1/3, suggesting a type of Sa or 
later. 

NGC 4586 This galaxy shows spiral arms and a normal 
bulge for an Sa galaxy. 
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Fig. 7. POSS I snapshots of the T=1H i-dehcient galaxies of Table 7. The box size is 6' and North points upwards. 



NGC 4606 This galaxy shows little traces of spiral arms 
and has an intermediate bulg e. NGC 4606 appears be- 
tween SO/a and Sa. However. iGavazzi et alJ |200Cfl fit 
a bulge to disk ratio of 1/20, suggesting a later type 
than Sa. 



Altogether, we suspect that 2 galaxies may be earlier than 
Sa: NGC 4064 and NGC 4424, both of which appear to 
have recently ingested smaller galaxies, and significant 
amounts of cold gas may have been heated up during the 
merger. 



Given that between types SO/a and Sa, the compilation 
of lBettoni et~ai1 l|2003h (left plots of their Fig. 2) indicates 
that the normalized HI content varies as logMgj™ ~ 
T/4 + est, an overestimate of 1 (0.5) in T will lead to an 
underestimate in H I-deficiency of 0.25 (0.125). The Hl- 
deficiency of NGC 4064 is sufficiently high that a correc- 
tion of 0.25 will still keep it among the 3 a cases using both 
DEFi and DEF 2 criteria. On the other hand, if the mor- 
phological type of NGC 4424 is SO/a (T = 0), it would 
still be Hi-deficient with criterion DEFi, but no longer 
among the 3 a cases with DEF2. 
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Note, however, that if most Hi-deficient outliers clas- 
sified as Sa's were in fact SO/a's, we would find ourselves 
with the new problem of justifying the presence of SO/a's 
so far away from the cluster core and apparently not 
within gro ups (see section 5.3). We als o have to keep in 
mind that iKoopmann fc Kennevl l)l998[) argued that spi- 
rals in the Virgo cluster could have also been misidentified 
as early-type objects, an opposite trend to the above sug- 
gested. 

4.5. Case by case analysis 

We now focus on each of the 13 galaxies, one by one. 

N GC 4064 The very low i nner rotation 
ijRubin. Waterman, fc Kennevl Il999l) suggests that it 
is a merger remnant (see NGC 4424 below), which 
may have been responsible for its H I-deficiency. 

NGC 4305 This galaxy, at 40 Mpc and with a large ra- 
dial velocity, is almost certainly a background galaxy. 
It has a major tidal companion (Table 8) and a blue 
color for its luminosity, assuming it is indeed far behind 
Virgo, confirming the tidal origin of its gas deficiency. 

NGC 4307 The origin of the H I-deficiency in this galaxy 
is difficult to ascertain, as it has 41% (48%) probability 
of lying within the sphere centered on the Virgo clus- 
ter with radius its 1 (2.5) times the virial radius. Its 
proximity with NGC 4356 suggests that it may have 
interacted in the past with that galaxy, and if it is 
indeed in the background of Virgo, as suggested by 
its Tully-Fisher distance, it would be blue for its lu- 
minosity, confirming the interaction hypothesis with 
NGC 4356. 

NGC 4312 This galaxy is very likely within the Virgo 
cluster. 

NGC 4356 The origin of the H I-deficiency in this galaxy 
is difficult to ascertain, as it has 19% (23%) probability 
of lying within the sphere centered on the Virgo clus- 
ter with radius its 1 (2.5) times the virial radius. Its 
proximity with NGC 4307 suggests that it may have 
interacted in the past with that galaxy, as confirmed 
by its blue color for its luminosity, assuming it lies in- 
deed in the background of Virgo. 

NGC 4424 The very low inn er rotation of this galaxy 
suggests a merger remnant llRubin et al.lll999]) . It al- 
most certainly lies in the foreground of the Virgo clus- 
ter. Moreover, as discussed in Sect. 14.41 NGC 4424 may 
be of type SO/a, in which case it would not be 3cr de- 
ficient in neutral Hydrogen. 

NGC 4438 This seemingly foreground galaxy is 
very probably a mem ber of the Virgo cluster. 
iKotanvi fc Ekersl l|l983h suggest that the HI- 
deficiency of this highly distorted galaxy is caused not 
by the tidal companion, which should increase the 
line-width, but by interaction with the intracluster 
medium of the Virgo cluster. Moreover, this galaxy 
has a very small radial exte nt of its HI gas, a s 
attested by the VLA map of ICavatte et alJ IjlflM) . 



Interestingly, at the distance of Virgo, this galaxy 
would be very blue, which is easily explained by the 
presence of its tidal companion. 

NGC 4492 This galaxy is probably in the background 
of the Virgo cluster, in which case it would be blue for 
its luminosity and morphological type. It may lie in a 
group that could include NGC 4307 and NGC 4356 and 
thus be suffering tidal encounters with these and/or 
other members of the group. 

VCC 1569 This faint galaxy probably lies in or very 
close to the Virgo cluster. 

M90 Among the 13 galaxies lis ted i n Tab le 7, this 
very blue Sab LINER (e.g. iKeell Il996|) deviates 
the most significantly from both the Tolman-Bondi 
locus (Fig. |2J and from the color-luminosity di- 
agra m (Fig. 13. Its distance h as be en controver- 
sial. IStauffer. Kennev. fc Younel (^586) suggest that 
the Tully-Fisher estimate is an underestimate be- 
cause the optical rotation curve is still rising, 
to th e limit of its measurement (see iRubin et alJ 
1999 ) and that similarly low distan ces derived by 
Cowlev. Crampton. fc McClurel (|l982j) using their line 
index / bulge luminosity relation may be erroneous be- 
cause the weakness of the nuclear lines is in contrast 
with the high mctallicity expected from the stellar pop- 
ulation of the bulge. 

M90 is extremely isolated in projected phase space. 
Among the galaxies with D < 13 Mpc, it appears 
associated with a group of ellipticals, none of which 
has accurate distance measurements, with sky coordi- 
nates associated with the central ellipticals with accu- 
rate distance measurements in the center of the Virgo 
cluster. Given its bright corrected total blue magni- 
tude of 9.63 (Table 7), it's absolute magnitude would 
correspond to Jl/ B = -20.24 at 9.4 Mpc or M B = 
—21.5 at 16.8 Mpc. The latter absolute magnitude is 
very bright for a spiral galaxy. This can be quanti- 
fied by comparing with the turnover luminosity, £*, 
of the field galaxy luminosity function. Unfortunately, 
the surveys used to derive the field galaxy luminos- 
ity function do not employ the Johnson B band used 
here. However, the 2 MASS survey has measured a to- 
tal if-band magnitude of K = 6.58 for M90. For a 
turnover absolute mag nitude l|Kochanek et alJ l200l|) 
M* K = -23.39 + 5 log h = -24.27 for h = 2/3 adopted 
here, the X-band luminosity of M90 corresponds to 
0.4 L* if M90 is at 9.4 Mpc but 1.3 L* if it lies at the 
distance of the Virgo cluster, which is a large lumi- 
nosity for a spiral galaxy. Moreover, at the distance of 
Virgo, M90 would be extremely blue for its luminosity 
(Fig. 

Finally, the contours of continuum radio emission seen 
from the NVSS survey taken at the VLA are com- 
pressed in the East-Northeast direction, which is op- 
posite to the direction towards the center of the Virgo 
cluster. This suggests that M90 may be bouncing out 
of the cluster, which would also explain the galaxy's 
negative radial velocity. 
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NGC 4580 This galaxy probably lies near but not in the 
Virgo cluster, and may possibly be bouncing out of the 
cluster, therefore losing its neutral gas by ram pressure 
stripping. It may be interacting with NGC 4586, and 
is somewhat blue for its luminosity, especially if it lies 
behind Virgo. 

NGC 4586 This galaxy probably lies near the Virgo 
cluster and may possibly be bouncing out of the clus- 
ter, therefore losing its neutral gas by ram pressure 
stripping. It may be interacting with NGC 4580, and 
is blue for its luminosity, especially if it lies within the 
Virgo cluster. 

NGC 4606 This galaxy is a very likely member of the 
Virgo cluster. It is likely to be tidally perturbed by its 
companion NGC 4607 (Table 8), as attested by its blue 
color for its luminosity assuming it indeed lies in Virgo, 
although it shows only modest signs of such tides. 



lully-i 1 isner distances may 

4. 6. Origin of the H l-deficiency in the outlying galaxies For exa mple, Effer et n ,lJ i 

The Hl-dcficicnt galaxies apparently lying well in front 
and behind of the Virgo cluster must have passed through 
the core of the cluster to have shed their gas by stripping 
by the ram pressure of the intracluster gas. As shown by 
iMamon et alJ l(2004h . galaxies cannot bounce out beyond 
1 to 2.5 virial radii, which according to their estimate of 
the virial radius of Virgo, amounts to 1.7 to 4.1 Mpc. 
Therefore, given a distance to the Virgo cluster of 16.8 
Mpc, we suggest that only galaxies with distances between 
12.7 and 20.9 Mpc could have passed through Virgo. 

One possibility, discussed in Sect. 14. II is that the true 
3D locations of the outlying H I-deficient galaxies are in 
fact within the Virgo cluster's virial radius or at most 
within 4 Mpc from its center, which is possible (although 
not likely) given the typical 20% relative distance errors 
for the spiral galaxies in Virgo, expected from the Tully- 
Fishcr relation, and from our comparison of the observed 
velocity distance relation with that derived from cosmo- 
logical TV-body simulations. Indeed, of the 13 H I-deficient 
spirals appearing in the outskirts of the Virgo cluster, 4 
to 5 (depending on how to place the distance to M90) 
are highly likely to be located in 3D within the virial ra- 
dius of the Virgo cluster (Table 7) and one other galaxy 
(NGC 4307) has one chance in two of being in the cluster 
or near it. Three others are at an angular radius > #i 00 
from M87, but appear to be within 2.5 virial radii of it 
in 3D. Finally, two galaxies (4356 and NGC 4492) in the 
direction of the cluster have low but non negligible proba- 
bility of being in or near the cluster, while the last two are 
definitely foreground (NGC 4424) and background (NGC 
4305) galaxies. Hence, between 4 and 11 of these 13 galax- 
ies may have passed through the core of the Virgo cluster 
and seen their neutral gas ram pressure swept by the hot 
intracluster gas. 

The presence of a few discrepant distance measure- 
ments suggests that the nearly gaussian parent distri- 
bution function for measurement errors in the distance 



modulus (or relative distance) has extended non-gaussian 
wings. 

It is also possible that highly inaccurate (at the 30% 
level or more) Tully-Fisher distance measurements may 
arise from physical biases. If galaxies arc H I-deficient be- 
cause of strong interactions with the cluster (through ram 
pressure stripping or tides) or with neighboring galax- 
ies (through collisional tides), their Tully-Fisher distances 
may be biased and/or inefficient. Indeed, our analysis of 
Sect. l4.3l indicates that, among the galaxies likely or possi- 
bly within the Virgo cluster, one (NGC 4438) shows clear 
tidal perturbations from a close major companion, while 
2 others (M90 and NGC 4606) are probably tidally per- 
turbed by close companions, and the one possible Virgo 
cluster member (NGC 4307) may have interacted with a 
fairly close major companion in the past. 

Moreover, there are general arguments, independent of 
these case-study tidal effects, that suggest that abnormal 
Tully-Fisher distances may be c aused by H I-deficiency. 



1986ft noted that H I-deficient 
spirals in the direction of the Virgo cluster tend to have 
small er line-widths for their if -band luminosity, and, sim- 
ilarly, iRubin. Hunter, fc Fordl l|l99l|) found that the max- 
imum rotation velocity of spiral galaxies within com- 
pact groups of galaxies tend to be smaller for their lu- 
minosity than for normal spirals. Such small line-widths 
or maximum rotation velocities imply that distances are 
underestimated. Using 2D (Fa bry-Perot) spectroscopy, 
iMendes de Oliveira et all <)2003|) have refined the analysis 
of IRubin et al.l and find that the lowest luminosity galax- 
ies in their sample (Mb = — 19.5±1) have lower maximum 
rotation velocities than field spirals of the same luminosity. 
Moreover, they find that compact group galaxies display 
a more scattered Tully-Fisher relation than that of field 
spirals. 

For galaxies lying in or very near the Virgo cluster, 
an alternative to ram pressure gas stripping is for tidally 
stripped gaseous reservoirs that prevent subsequent gas in- 
fall to the disk, which quenches sub sequent star formation 
1 Larson. Tinslev. fc Caldwelll ll980). However, if a spiral 
galaxy lies in a cluster, it may well be on its first infall, 
provided that the effective ram pressure stripping of its 
neutral gas quickly converts it into an SO galaxy. It seems 
reasonable to expect that tidally stripped gas reservoirs 
will not have time to have a large effect on the gas content, 
so that such galaxies with tidally stripped gas reservoirs 
will not be as severely H I-deficient as the 13 3 a cases of 
Table 7. 

The H I-deficient galaxies that are truly outside of the 
cluster may have lost their gas by the tides generated dur- 
ing close encounters with other galaxies within a group 
or small cluster of galaxies. Indeed, the efficiency of tidal 
perturbations of galaxies from close encounters with other 
galaxies is stronger in gr oups and sma ll clusters than in 
rich clusters (see Fig. 5 of lMamonl2000|) . basically because 
the rapid motions in rich clusters decrease the efficiency of 
tides. Three of the four H I-deficient galaxies at 28 Mpc lie 
within a significant concentration of galaxies (Figs. ^ and 
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0J) that ought to create such collisional tides, even with 
the estimates of distance using the more precise radial ve- 
locities (for these distant objects, see eq. |T2]1. 

The 3 galaxies outside the virial angular radius (NGC 
4064, NGC4580 and NGC4586) all probably lie closer than 
2.5 virial radii from M87, and may be in the process of 
bouncing out of the cluster after having lost their H I gas 
through ram pressure stripping by the intracluster gas, if 
the maximum rebound radius is closer to the upper limit of 
2.5 virial radii found bv lMamon et al , than of their 

lower limit of 1 virial radius. One of these galaxies (NGC 
4064) is a candidate for a recent merger as witnessed by 
its very low inner rotation, while the other two (NGC 
4580 and NGC 4586) may have interacted in the past, 
given their proximity in phase space. In the case of recent 
merging, the neutral gas may have been shock heated with 
that of the galaxy being swallowed up, and should shine 
in Ha or X-rays depending on its temperature. 

Of the 4 galaxies likely or certainly far outside the 
Virgo cluster (in 3D), one (NGC 4305) may be tidally per- 
turbed by a close major companion, one (NGC 4424) has 
suffered a recent merger (given its very low inner rotation 
velocity), one may have interacted with another galaxy 
(NGC 4307) in the past, and one (NGC 4492) may be a 
misclassified SO/a (with T = 0), implying a lower expected 
HI content, hence a less than 3cr H I-deficiency upon the 
estimator DEFi and just at 3 a upon DEF? , alth ough the 
bulge/disk ratio derived bv lGavazzi et al.1 <)2000h suggests 
that the galaxy is indeed an Sa. 

Although we can provide at least one explanation for 
each case of Hl-dcficicncy in the 13 spiral galaxies that 
1) appear to be located in the foreground or background 
of the Virgo cluster or 2) lie just outside the projected 
virial radius of the cluster, there is certainly room for an 
improved analysis. In particular, accurate distances ob- 
tained with Cepheids would be most beneficial, as well as 
VLA maps for 11 of the 13 galaxies to distinguish tidal 
stripping from ram pressure stripping. 
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Appendix A: Median and half-width of the 

likelihood function for a set of measurements 
with different errors 

For a single measurement of variable X of value x with 
uncertainty cr^, assuming gaussian errors, the probability 
of measuring Xi is 

\2- 



Pi(xi\x) = 



1 



2tt<j, 



exp 



2aj 



(A.1) 



The likelihood of measuring a set of {xi} given a true value 
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From Bayes' theorem, the probability of having a true 
value x is 



P(x\{ Xi }) 



f(x) 



£ ({xi}\x) 



(A.4) 



where f(x) is the a priori probability of having a true 
distance x and g({xi\) is the a priori probability of mea- 
suring the set of {x^. If the measurements are made with 
a selection function near unity, then g ~ 1. 

For simplicity, one can assume a uniform f{x). The 
probability of the true X being less than some value x is 
then 



P{< x) 



Hoo^ijxiWx) dx 

^ exp [-x 2 /2 gj Vj +xJ2 t Xi/crf] dx 
exp [-a?/2 J2i y°l + x Ei xjaf] dx 
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where 



Ei i/of 



(A.6) 



Now erfc(zVV2)/2 = 0.16,0.5,0.84 for v = -1,0, +1, so 
that the median X is 

T _ E, x 1 /aj 

50 " E7W ' ( } 

while the width equivalent to the standard deviation of a 
single gaussian is 
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Since C is an exponential function of x (eq. [A.2| )„ the 
mean and mode will be equal to x 50 . 

Alternatively, one could adopt an underlying physical 
model for fix). If X is the distance modulus, for which the 
gaussian p.d.f. of equation i|A.l|) should be adequate, one 
could use a Virgocentric model, where the probability f[x) 
stems from a density distribution around Virgo, using for 
example an NFW model. If the angle between the galaxy 
and Virgo, as seen by the observer, is 9, one has 

r 2 = D 2 10 {dcx(-0.4d) + dcx(-0.4d y ) 

-2 cos(9dex[-0.2(d + dy)]} , (A.9) 

where d and dy are the distance moduli to the galaxy and 
to Virgo, respectively, and Z?io = 10 pc. Hence, one has 

f id) K p ( r) *L } (A.10) 

dr 1 dr 2 
dd ~ 2r ~dd 

0.2 lnl0£> 2 

r 

x [cos6>dex(-0.2(d + d y )) - dex(-0.4d)] .(A.ll) 

The problem with this approach is that the underly- 
ing model, i.e. p(r), is uncertain: Should one take a single 
NFW model at the location of M87? Should one make 
p(r) uniform at a radius where PNFw( r ) reaches the mean 
mass density of the Universe? Or better, should one re- 
place this uniform background with a collection of other 
halos, as expected in the hierarchical scenario? Moreover, 
because the model requires a distance to M87, one would 
need to iterate using a first guess on the distance to M87 
(perhaps based upon the estimate with a uniform f(x)). 
These questions led us to assume in this paper a uniform 
distribution for f(x). 
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